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The thermal evolutions of the Moon, Mars, Venus, and Mercury are calculated theoretically 
starting from cosmochemical condensation models. An assortment of geological, geo- 
chemical, and geophysical data are used to constrain both the present-day temperatures 
and the thermal histories of the planets’ interiors. Such data imply that the planets were 
heated during or shortly after formation and that all the terrestrial planets started their 
differentiations early in their history. 

Initial temperatures and core formation play the most important roles in the early 
differentiation. The si 2 e of the planet is the primary factor in determining its present-day 
thermal state. The Moon, smallest in size, is characterized as a differentiated body with a 
crust, a thick solid mantle, and an interior region which may be partially molten. It is 
presently cooling rapidly and is relatively inactive tectonically. 

Mercury, which probably has a large core, may have a 500-km-thick solid lithosphere 
and a partially molten core, if it is assumed that some heat sources exist in the core. If 
this is not the case, the planet's interior temperatures are everywhere below the melting 
curve for iron. The thermal evolution is dominated by the core separation and the high 
conductivity of iron that makes up the bulk of Mercury. 

Mars, intermediate in size, is assumed to have differentiated an Fe-FeS core. While the 
formation of an early crust is evident, large-scale melting and differentiation of the mantle 
silicates has occurred at least up until 1 b.y. ago. Present-day temperature profiles indicate 
moderate tectonic activity at the present time. 

Venus is characterized as a planet not unlike the Earth in many respects. Core formation 
has occurred probably during the first billion years after the formation. The present-day 
temperatures indicate a partially molten upper mantle overlain by a 100-km-thick litho- 
sphere and a molten Fe-Ni core. We can expect that today Venus may have tectonic 
processes similar to the Earth’s. 


The evolution of a planetary body is very 
strongly controlled by its thermal regime. 
Differentiation, volcanism, tectonic activity, 
and even magnetic field are all processes 
controlled by the internal temperatures and 
the thermal history. The calculation of the 
thermal evolution, however, requires strong 
constraints if it is to be realistic. In recent 
years significant new data from lunar and 
planetary missions have become available. 
For the Moon these have provided constraints 
both on the present-day temperatures and 
the past history of the thermal state of the 


lunar interior. For Mars, Venus, and Mer- 
cury, the thermal constraints are primarily 
qualitative. The purpose of this paper is to 
review these data and the thermal evolution 
models of the Moon and the terrestrial plan- 
ets. 

An extensive number of studies on thermal 
history of the Moon and planets have been 
carried out (refs. 1 through 34). 

In the early studies, the Moon and the 
planets were generally modeled with chon- 
dritic compositions and radioactivity. With 
the successful lunar missions and new geo- 
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chemical data, it became necessary to revise 
the constraints used in lunar thermal calcu- 
lations. The models generated since 1971 
have incorporated these new data and in 
general give similar results. Generally, cal- 
culations for the other planets reflect this 
new knowledge about the chemistry of the 
solar system. 

In this paper, then, we present thermal 
evolution models for the Moon, Mars, Venus, 
and Mercury for a comparative study. We 
use similar computational techniques and 
thermal factors associated with melting, dif- 
ferentiation, and core formation. These are 
discussed first. The thermal evolution of the 
Moon, with appropriate constraints, is then 
presented, followed by a discussion of the 
thermal evolution of the planets. 

Computational Procedure and 
Model Parameters 

For the theoretical calculations of thermal 
evolution in a planetary body, it is necessary 
to specify the initial conditions, heat sources, 
heat transport mechanisms, and a number 
of other parameters such as specific heat, 
density, and melting temperatures, as well 
as variations of these with temperature and 
pressure inside the planet. Generally, not all 
of these parameters can be determined and 
theoretical calculations cannot define thermal 
models uniquely. 

The first input in theoretical calculations 
is the initial conditions assumed. Included 
in this is the heating due to accretion, tidal 
forces, induced electric currents, and short- 
lived radioactivity, although the relative im- 
portance of each may vary. Next is the 
definition of the long term heat sources that 
are important throughout the age of the 
planet. These include heating due to the decay 
of radioactive isotopes and the possibility of 
long term heating due to tidal dissipation 
and electrical currents. Finally, using the 
above inputs combined with other physical 
parameters, the conservation of energy equa- 
tion must be solved numerically with the 
addition of certain perturbations designed 


to simulate melting, differentiation, and con- 
vection. 

In this section we first describe the input 
parameters and then discuss the computa- 
tional procedures involved in solving the 
thermal evolution problem. 

MODEL PARAMETERS 
Initial Temperature 

The initial temperature of a planet is a 
function of both its mode of formation and 
its immediate environment during and 
shortly after its origin. As will be discussed 
in the following sections, constraints on lu- 
nar and planetary thermal evolution favor 
initially hot and possibly molten models. 
This is especially true for the Moon. Energy 
sources which might have provided the heat 
to raise the temperature of the planets to 
solidus or near-solidus levels include gravi- 
tational energy due to accretion, tidal 
dissipation, solar wind flux, short-lived ra- 
dioactivity, and heating due to adiabatic 
compression. The temperature rise due to 
adiabatic compression amounts to only a few 
tens of degrees in the Moon, but increases 
with the increasing size of the planet. The 
effect of short-lived radioactivity depends on 
the time between nucleogenesis and lunar 
formation. Also, Schramm et al. (ref. 35) 
have shown that AF«, the most likely heat 
source in the early history of the solar system 
(ref. 36), had little effect on the initial heat- 
ing of the Moon and meteorites. 

Tidal dissipation and solar wind flux could 
have been effective under special circum- 
stances. The amount and the distribution of 
heat is difficult to quantify in both cases. For 
the Moon, as an example, the effects of tidal 
dissipation would be small (ref. 37), unless 
the Moon’s orbit was originally much closer 
to the Earth. This might have occurred 
shortly after a capture event or after the 
Moon formed in close proximity to the Earth 
(ref. 38). Hallam (ref. 39) has calculated 
that tidal dissipation, even with Q = 10, 
raises the internal temperature by only 50- 
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lOOK in 1.5 b.y. Even if the Moon was once 
close to the Roche limit, it would recede so 
rapidly that the effect of tidal heating would 
be negligible. 

Heating by electric currents produced by 
a unipolar generator driven by the solar wind 
might have been significant if the Sun had 
passed through a T-Tauri stage (ref. 40). 
This effect is difficult to calculate and include 
explicitly in any thermal model. However, 
combinations of all the above effects may 
contribute to the initial temperatures. We 
considered these effects by means of a “base” 
temperature added to accretional thermal 
profiles discussed below. 

The most promising initial heat source for 
the terrestrial planets is the gravitational 
energy of accretion. Classically this effect is 
described by the equilibrium of added gravi- 
tational potential energy with black body 
radiation and internal heat; 

= ea(T(r)*-Ti)dt 
'' +pCAT(r)-T,)dr (1) 

where t is time, p is the density of the ac- 
creting particles, G is the gravitational con- 
stant, M (r) is the mass within radius r, (7 is 
the Stefan-Boltzmann constant, Cp is the 
specific heat, and e is the emissivity (taken 
to equal 1). T(r) is the temperature at ra- 
dius r, and Tt, is the base temperature or 
the temperature of the accreting particles, 
including any effects of short-lived radioac- 
tivity. By specifying the accretion rate dr/ 
dt and the physical parameters p, Cp, and 
Tb, temperature versus depth for an accreted 
body can be calculated from equation (1). 
An accretion rate suggested by Hanks and 
Anderson (ref. 15) and Mizutani et al. (ref. 
41) is used in this paper and is given by: 

= ct^ sin y t (2) 

The constants c and y are determined from 
the planet’s final radius and the duration of 
accretion. Although the accretion rate is 
difficult to establish, it is a fair statement of 
the presumption that the rate grows as the 
accreted body grows, and that the rate then 


tapers off as primary material is depleted. 

The accretion time is considered arbitrary, 
except that it is assumed to be longer for 
larger planets. Generally the accretion time 
and Tb are adjusted to provide an initial 
temperature profile consistent with the con- 
straints that require initial differentiation. 
For the Moon, it can be shown from equa- 
tions (1) and (2) that the temperature will 
never reach the basalt solidus for accretion 
times greater than 1000 years, unless the 
emissivity of the lunar surface is decreased. 
This conclusion was reached independently 
by Mizutani et al. (ref. 41) from a more de- 
tailed study of the accretion process. 

Heat Sources 

The most important long term heat source 
in thermal history calculations is the heat 
generated from the decay of the long half- 
life radioactive isotopes U*®®, Th“®, 

and K^®. Abundances of these isotopes have 
been measured in many Earth rocks and re- 
turned lunar samples (refs. 42 through 48, 
and others), and indirectly by lunar orbital 
y-ray spectrometers (ref. 49). While surface 
concentrations vary greatly, we may conclude 
that these values must be much higher than 
those in the interior of the Earth or the Moon 
(ref. 50, and others). The distribution of the 
isotopes in space and time are difficult to 
constrain, however, due to complex melting 
and crystal-liquid fractionation processes. 

Figure 1 shows total K and U abundances 
measured for Apollo samples (including crys- 
talline rocks, soils, and breccias), some me- 
teorites, and terrestrial rocks. Figure 2 shows 
Th/U versus U for the same suite of mea- 
surements. Although there are small and sys- 
tematic deviations, the K/U ratio for lunar 
material averages around 2000, compared 
with 10 000 for terrestrial rocks and 80 000 
for chondrites. The low lunar K/U ratio is 
consistent with general depletion of volatile 
elements (refs. 42 and 52). From this, it is 
clear that neither terrestrial nor chondritic 
abundances can be used as the bulk radio- 
activity for the Moon. Furthermore, using 
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Uranium, ppm 

Figure 1 . — Potassium and uranium concentrations in selected meteorites, lunar samples, and terrestrial rocks. 
Lunar rocks are separated into three categories: crystalline rocks represented by filled symbols, breccias 
shown by half -filled symbols, and fines shown by empty symbols. Points for average choridrites, howardites, 
and eucrites are from Mason (ref. 51 ). Sources for lunar values are given in the text. 


Rb abundances for a model of the primitive 
lunar material (ref. 18), we can also reject 
chondrites. The Th/U ratio is consistently 
between 3.0 and 4.0 for terrestrial, lunar, and 
meteoritic material. 

It is essential, then, to determine or esti- 
mate the bulk radioactive concentrations for 
each of the planets. In the case of the Moon, 
given that chondrites and terrestrial rocks 
are not good approximations of the bulk ra- 
dioactivity, we turn to achondrites — specifi- 
cally howardites and eucrites. A genetic 


relationship between achondrites and the 
Moon has been suggested by several investi- 
gators (refs. 53, 54, and 55). Howardites 
have, on the average, a uranium concentra- 
tion of 35 ppb, with K/U = 9800, while brec- 
ciated eucrites, enriched in refractory 
elements, have U = 130 ppb, with K/U = 
2300 (ref. 51). It is interesting to note that 
the K/U ratio for eucrites is close to that of 
the lunar material, while the ratio for 
howardites is somewhat higher. 

If we take K/U for the bulk radioactive 
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Uranium , ppm 


Figure 2. — Th/U ratio versus uranium concentration for selected lunar rocks and meteorites. Lunar rocks 
are divided into crystalline rocks, breccias, and fines, with the same notation as in figure 5. Points for 
average chondrites, howardites, and eucrites are from Mason (ref. 51). Sources for lunar values are 
given in the text. 


concentration of the Moon to be 2000 and 
Th/U to be 3.5, then the average radioactiv- 
ity of the Moon is a function of the bulk 
uranium content. Given the measured Apollo 
15 and 17 heat flows of about 30 ergs/cm^-s 
(refs. 56 and 57), Toksoz and Solomon (ref. 
30) have shown that the bulk uranium con- 
centration must be 60 ± 15 ppb to match the 
heat flow value. In this paper, then, we let 
U = 60 ppb, with K/U = 2000 and Th/U 
= 3.5 for the Moon. 

For Mars, Venus, and Mercury, we will 
utilize values based on condensation models 
of the solar nebula (refs. 54, 58, and 59). 
These are shown in table 1. In the calcula- 
tions we consider variations for individual 
cases as discussed in later sections. However, 
we expect that for Venus terrestrial abun- 
dances are a good approximation. For Mars 
we use values between terrestrial and chon- 
drite abundances (ref. 32). For Mercury the 
uranium values are close to those of the Moon, 


but the planet is most likely lacking in 
(ref. 33). The relevant decay constants, de- 
cay energies, and isotopic abundances are 
taken from Clark (ref. 60). 


Distribution of Heat Sources 

In most cases, the thermal models pre- 
sented in this paper assume an initially 
homogeneous distribution of heat sources. 
However, the constraints on the internal tem- 
peratures rule out, at least for the Moon, a 
present-day uniform distribution of heat 
sources. When melting occurs, one would ex- 
pect the magma to be enriched in U, Th, and 
K and to subsequently transfer these heat 
sources to the surface. To account for this, 
at discrete time steps the heat sources from 
the molten zones are differentiated upwards, 
leaving primordial radioactive abundances 
in any solid region below the zone. Heat 
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Table 1 . — Parameters Used in Thermal History Calcidations 



Mercury 

Venus 

Moon 

Mars 

Radius (km) 

2439 

6050 

1738 

3389 

Mean Density (g/cm“) 

5.43 

5.25 

3.34 

3.96 

C/Ma= 

— 

— 

0.394 

0.377 

Surface Temperature (°C) 

77 

327 

-20 

-40 

Uranium, ppb 

44 

31 

60 

3 

CO 

K/U ratio 

(1> 

10,000 

2000 

10,000 “> 

Th/U ratio 

4 

4 

3.6 

4 

Specific Heat (silicates) 
(J/g:°C) 

1.2 

1.2 

1.2 

1.2 

Heat of Fusion (silicates) 

(J/g) 

400 

400 

400 

400 


Note: (1) See text. 


sources differentiated toward the surface are 
allowed to decrease exponentially with depth 
so that 

H{r) =Aoexp|^— (3) 

where .4o is the surface abundance, R is the 
radius, and H{r) is the radioactive abun- 
dance at radius r. Such a distribution has 
been formulated for the Earth (ref. 61). The 
skin depth h decreases with time and is ad- 
justed so that the present-day surface ra- 
dioactivity, if known, is approximated. 

Thermal Conductivity and 
Other Parameters 

In a differentiating planet, heat is trans- 
ported by conduction and by convection. At 
temperatures far below solidus, the role of 
convective flow is small and heat transport 
is governed by thermal conductivity. This 
quantity is temperature-dependent and is 
taken to be the sum of two terms, lattice (or 
phonon) conductivity and radiative conduc- 
tivity. The relative dominance of each term 
is of crucial importance in solving the con- 
duction problem of the planets. Although 
effective thermal conductivities have been 
measured for terrestrial and lunar materials 
(fig. 3), these do not always represent the 
bulk properties of the Moon or planetarv in- 


teriors. MacDonald's (ref. 8) theoretical 
formulation of conductivity as a function of 
temperature assumed the opacity to be con- 
stant, resulting in a cubic temperature de- 
pendence of the radiative term and, thus, its 
dominance over all temperatures. However, 
Schatz and Simmons (ref. 64) proposed that 
for a hypothetical dunite, radiative transfer 
is suppressed by grain-boundary extinction, 
resulting in an increase in the opacity with 
temperature. At low temperatures, then, the 
effective conductivity is dominated by the 
lattice term that decreases with increasing 
temperature. At higher temperatures, the 
radiative term, increasing roughly linearly 
with temperature, becomes dominant. This 
model has been used in previous thermal his- 
tory calculations at MIT (refs. 23 and 30). 
For the Moon, however, upon examination of 
figure 3 and values of conductivity in many 
previous lunar models, one could just as well 
take the effective conductivity to be a con- 
stant. In some of the lunar models and for 
Mercury’s mantle we shall treat it thus, with 
K = 0.45 X 10® ergs/cm-s °C, while other cal- 
culations assume the model of Schatz and 
Simmons (ref. 64). The total conductivity, 
K, is given by 

where = lattice conductivity and Kr — 
radiative conductivity. As functions of tem- 
perature, T, these are given by 
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T. °C 



Figure 3 . — Thermal conductivity of selected mate- 
rials. Sources of measured data (solid lines) are 
Cremers (ref. 62) for Apollo 11 fines, sample den- 
sity 1.6A g/cm”; Horai et al. (ref. 63) for Apollo 
11 basalt; Murase and McBirney (1970) for syn- 
thetic Apollo 11 basalt; Schatz and Simmons (ref. 
6i) for sintered, poly crystalline forsterite and for 
single crystal olivine (Foss Fan) cmd enstatite. 
Theoretical (dashed) curves include those of Mac- 
Donald (ref. 8) for two different values of the 
opacity e, and one proposed by Schatz and Simmons 
(ref. 62) for poly crystalline olivine of approximate 
composition FouFatt. The Schatz and Simmons 
curve, labeled S above and K = constant, was em- 
ployed in this paper. 


^ __ 4.184 X 10'^ 

^ 30.6 + 0.21T 

and 

^ _ (0 for r < 500°K 

~ |230(r-500) for T > 500°K 

Other physical parameters required for 
the calculations are the melting curve, sur- 
face temperature, specific heat, heat of 
fusion, and density. A summary of the 
parameters common to all models is given in 
table 1. Appropriate melting curves are used 
for each planet, taking into account the com- 
position and pressure. 

COMPUTATIONAL TECHNIQUES 

In this section we describe the thermal 
calculations for a spherically symmetric 
planet, taking into account conduction, melt- 


ing, simulated convection of molten material, 
and differentiation. 

Conductive Models 

Assuming that conduction is the only mode 
of heat transfer in solid material (an approx- 
imation to convection will be discussed later) , 
thermal evolution models are calculated us- 
ing the finite difference solution of the heat 
conduction equation: 

if) + w 

where Cp is the specific heat, T is the tem- 
perature, r is the radius, K is the thermal 
conductivity, and H (r,t) is the heat source 
term. The latter two quantities were dis- 
cussed above. 

The details of the calculations are de- 
scribed fully in the appendix of Toksoz et al. 
(ref. 23) and are stated here briefly. The 
finite difference analog of equation (4) 
which conserves heat flux is 



{K H 

where p = Ar, r — rvp, t = mAt, « = At/pCp, 
and T w= T(np, mAt). The stability condition 
relating the time increment and the grid spac- 
ing is 


(A„:,+ A„") 


(T„”-n) 


n 


At < 


CppAr^ V? 



( 6 ) 


where Amar is the maximum value of the 
conductivity at a given time step. At is com- 
puted at each iteration using equation (6), 
with a factor of 4 instead of 2 in the de- 
nominator. A grid spacing, Ar, of 20 km is 
used. 

Boundary conditions for the problem re- 
quire the surface temperature of the planet 
to be constant, and T” - ■ T”^, where T”^ is the 
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temperature at ri = Ar and T" is the tem- 
perature at r2 = ri + Ar. The solution starts 
at r — Ar to avoid a singularity at r = 0. 
Surface heat flux is calculated by adding the 
heat flow at the boundary between the upper 
two finite difference cells to the steady-state 
heat flow due to heat production in the upper- 
most half cell. 


Simulation of Melting and Convection 

The effects of melting and convection of 
molten material on temperature within the 
body can be modeled following the technique 
described by Reynolds et al. (ref. 65). In 
this method the temperature at each grid 
point at each time step is compared to the 
melting temperature at the respective grid 
point. If the temperature is greater than the 
melting point, the difference divided by the 
density and specific heat gives the heat 
equivalent, If aH is less than or equal 
to the heat of fusion, AHf, the material at 
that grid point is considered partially molten, 
and the temperature is set equal to the 
melting point. This process is continued at 
each iteration until the accumulated AH is 
greater than AHf. At this point in time, the 
material is considered to be completely 
molten and the temperature is allowed to rise 
above the solidus. The same process allows 
completely molten material to solidify, re- 
leasing heat equal to AHf. 

In simulating convection of molten ma- 
terial, the same technique as outlined above 
is followed, with the exception that the tem- 
perature of the completely molten material 
is constrained to the melting curve. Any in- 
crease in temperature that would raise the 
temperature above the melting point is con- 
verted to its heat equivalent, transferred 
upward to the next grid point, and consid- 
ered to be a temperature rise due to con- 
vection. The conversion to heat equivalents 
is necessary due to the increase in the vol- 
ume element with increasing radius. The ad- 
vantage of this approximation to convection 
is that it does not require the inclusion of 
mass transfer terms in the heat conduction 


equation. This scheme is not limited to con- 
vection of completely molten material. Con- 
vection of partially molten material is easily 
approximated by reducing the heat of fu- 
sion. Convection by solid-state creep requires 
techniques different from those described 
above. Various attempts at approximating 
solid-state convection and determining its 
relative effect on the thermal evolution of 
the Moon and planets are discussed in the 
following section. 


SOLID STATE CONVECTION 

The question of whether or not thermal 
models calculated assuming conduction mech- 
anisms are unstable to solid-state convec- 
tion motions has been discussed in many 
papers, most recently by Tozer (refs. 24 and 
25), Turcotte et al. (ref. 66), Cassen and 
Reynolds (refs. 29 and 67), and Schubert 
et al. (ref. 34). Much of this work has been 
motivated by Schubert et al. (ref. 68), Us- 
ing a number of conduction temperature pro- 
files for the planets, these investigators 
applied a linear stability analysis and con- 
cluded that thermal convection would occur 
in each case. Earlier, Runcorn (refs. 69 and 
70) attempted to explain departures of the 
lunar surface from hydrostatic equilibrium 
by employing convection. Whether or not 
convection has played a dominant role in the 
evolution of the planets, however, is still de- 
batable. The question is not easily answered, 
since the presence (or absence) of convec- 
tion is a complex function of the temperature 
and the rheology of the planetary material. 
Futhermore, in the case of the Moon, the 
characteristics of the mare basins and their 
possible filling with basalts by gravitational 
upwelling of magma implies the presence of 
more complicated geometries than the sim- 
ple cellular convection pattern postulated by 
some investigators. Such upwelling type 
heat transport is well approximated in ou: 
calculations described earlier. Also, the at- 
tenuation of shear waves and the presence 
of deep moonquakes in the lunar interior, as 
discussed in the following section, present 
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difficulties for those who hypothesize that 
solid-state convection occurs in the Moon. 

The strongest advocate for the dominance 
of solid-state creep in the Moon and planets 
has been Tozer (refs. 24, 25, and 26). In all 
his papers, Tozer has postulated that con- 
vection is at least as important as conduc- 
tion for any planetary body of radius greater 
than a few hundred kilometers. This theory, 
based on the assumption that at any finite 
temperature an arbitrarily small nonhydro- 
static stress produces some permanent defor- 
mation of the material, suggests that at 
some viscosity, or range of viscosities, this 
motion, or rate of creep, becomes “fast” 
enough to provide a very efficient method 
of heat transfer. 

In this paper we do not go into detailed 
convection calculations. In some cases, how- 
ever, we include models that approximate 
solid-state creep by arbitrarily raising the 
thermal conductivity of material with tem- 
peratures above a certain threshold value. 

Thermal Evolution of the Moon 

CONSTRAINTS ON 
THERMAL EVOLUTION 

Thermal evolution calculations generally 
require boundary conditions that would in- 
dicate directly or indirectly the temperatures 
at certain times inside the planetary body. 
In addition to these, data are needed to 
specify material behavior, heat transport 
mechanisms, and heat sources. Normally, in- 
sufficient data and constraints are available. 
Thus, calculations generally are carried out 
with a set of initial conditions and assump- 
tions, and the conditions are modified if the 
temperature models fail to satisfy available 
constraints. 

Data and samples from Apollo missions 
imposed some strong constraints on the ther- 
mal evolution of the Moon. Some of these, 
such as heat flow, electrical conductivity, 
seismic velocities, attenuation, and tecto- 
nism, pertain to present-day temperatures. 
Some others, such as the presence of a dif- 


ferentiated lunar crust and chronology of lu- 
nar igneous and magmatic activity, indicate 
the thermal state early in the lunar history. 

Lunar Heat Flow 

Two heat flow measurements have been 
made on the Moon, one each at the Apollo 15 
(Hadley Rille) and Apollo 17 (Taurus- 
Littrow) sites. The experiments and their 
results are described by Langseth et al. (ref. 
57). The heat fluxes have been determined 
by measuring temperature gradients and 
thermal conductivities in bore holes to a 
maximum depth of 234 cm. In that depth 
range, temperature gradients in the lunar 
regolith are between 1.3 and 1.7 °K/m, and 
the thermal conductivity is in the range 
1.7 X 10-^ and 2.0 X 10-^ W/cm-°K. The 
heat flow through the surface based on the 
measured temperature gradients and con- 
ductivity is 31 ergs/cm*-s at the Apollo 15 
site and 28 ergs/cm^-s at the Apollo 17 site 
(ref. 57). In each case the estimated un- 
certainty is about ± 20 percent. 

The Taurus-Littrow and Hadley Rille sites 
are located in the basalt-flooded valleys at 
the rims of the Serenitatis and Imbrium ba- 
sins, respectively. There are no direct data 
to suggest that the thermal characteristics 
of these sites are any different from other 
parts of the Moon. Thus, an average heat 
flow value of about 30 ergs/cm^-s will be 
adopted for the Moon in these calculations. 
This value is about one-half the average heat 
flow of the Earth (63 ergs/cm^-s). If we 
assume steady-state conditions for both bod- 
ies and take into account the differences in 
the radii (average radii: Rifoon'= 1738 km, 
REarth = 6371 km), surface area (Smood/ 
<S*Earth “ 0.074), and volume (FMooa/TEarti. 
= 0.02), we find that the lunar heat flow is 
higher by about a factor of 2 compared to 
the Earth. As was discussed earlier in this 
paper, this implies higher heat generation 
and radioactive sources inside the Moon than 
inside the Earth. 

Another important measurement arising 
from the heat flow experiments is the value 
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of the near-surface temperature. At a depth 
of about 70 cm, the measured average am- 
bient lunar temperature is 253 K. This will 
be used as the surface temperature in our 
calculations. 

Electrical Conductivity and 
Inferred Temperatures 

Three magnetometers deployed on the lu- 
nar surface, one each at the Apollo 12, 15, 
and 16 sites, have recorded the magnetic 
fields induced in the Moon by large-scale 
transient events. In addition to these, the 
ambient and time-dependent magnetic fields 
in the lunar environment are measured 
by the Explorer 35 satellite magnetometer 
which is orbiting the Moon. 

The analysis of the transient fields re- 
corded simultaneously by the Explorer 35 
magnetometer and the surface magnetome- 
ter has made it possible to calculate the elec- 
trical conductivity inside the Moon (refs. 71 
through 74). The Moon responds to tran- 
sient magnetic events which induce eddy 
currents in the lunar interior. The response 
or transfer function is computed from the 
induced field at the Moon’s surface and the 
forcing (source) function recorded by 
the Explorer 35 magnetometer. 

These results are interpreted by calculat- 
ing the theoretical response of the Moon 
with radially varying conductivity. Al- 
though there are some limitations of the 
theoretical calculations and some difficulty 
in matching all the available data, the re- 
sults definitely show an increase of electrical 
conductivity with depth into the Moon. To 
a depth of about 1000 km the conductivity 
profile can be characterized by three layers 
(ref. 72) or by a greater variety of models 
including current layers (ref. 74). The re- 
solving power of the data is not sufficient to 
detail the conductivity below this depth. 

Conductivity depth profiles given by Dyal 
et al. (ref. 72) and the three-layer model 
(a good representation for the average mod- 
els) Tf Sonett et al. (ref. 74) are shown in 
figure 4. Basically, the models require very 


low conductivity (cri < 10 ® mho/m) for the 
outer 50 km of the Moon. For the interme- 
diate layer (1100 < R < 1700 km; B is ra- 
dius), the conductivity is in the range of 
a 2 = 1 to 7 X 10-^ mho/m. For the deeper in- 
terior there is a greater discrepancy. While 
Sonett et al. (ref. 74) give a value of about 
ITS ^ 10-® mho/m for R < 1200 km from the 
sunlit side data, Dyal et al. (ref. 72) give 

~ 10-^ for R < 1000 km. 

From the conductivity profiles, the tem- 
peratures inside the Moon can be estimated 
if the composition and temperature depen- 
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Fi^re 4 . — Lunar electrical conductivity profiles of 
Sonett et al (ref. 7i) and Dyal et al. (ref. 72) and 
some estimates of present-day temperatures in the 
Moon inferred from these measurements. Included 
are the temperature curve of Sonett et al. (ref. 7S), 
an interpretation of the Sonett et al. ( ref. 73) con- 
ductivity distribution using the temperature- 
conductivity relation of an olivine with little or no 
Fe‘* (ref. 28), and the estimates by Dyal et al. 
(ref. 71 ) based on their conductivity models. 
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dence of conductivity are known. For a 
given composition, conductivity, ar, increases 
with temperature : ar = ao exp ( — A ) , where 

kT 

A = activation energy, fc = Boltzmann con- 
stant, T — absolute temperature, and ao is a 
constant dependent on the material. If 
olivine were assumed for a composition and 
a lower bound on the conductivity of olivine 
is used, a temperature of about 1000 °C is 
estimated for the Moon at about 700 km 
depth (refs. 72 and 74). If newer laboratory 
data on pyroxenes and olivines — the most 
likely candidates for mantle composition 
(refs. 28, 75, and 76) — are used, higher tem- 
peratures (i.e., 1000°C at 400 km and 1400° 
at 800 km depth) are estimated (refs. 28 
and 77) . Some bounds on temperature models 
are shown in figure 4; in all cases tempera- 
tures are below the solidus to a depth of 
about 600 to 700 km. At greater depths, 
temperatures either approach or exceed the 
solidus, or some changes in the composition 
increase the conductivity. Another feature 
of this deep region is that both conductivity 
and temperature profiles flatten. 

Structure of the Lunar Interior and 
Seismic Wave Attenuation 

Seismic studies of the Moon have re- 
vealed much valuable information on the 
structure and the rheological properties of 
the lunar interior. The presence of a layered 
structure indicates differentiation and melt- 
ing in the past. The attenuation characteris- 
tics of seismic waves may be explained in 
terms of increased temperatures or partial 
melting at the present. 

The seismic velocity structure for the 
outer 150 km of the Moon has been deter- 
mined from the analysis of manmade im- 
pact seismograms recorded by the Lunar 
Seismic Network (refs. 78, 79, and 80). This 
structure shows the presence of a lunar 
crust overlying a mantle. Over the mare the 
crust may be two layered (the upper layer 
being mare basalt) , with a total thickness of 
about 60 km. In the highlands and on the 
backside, the crust may be somewhat thicker. 


The compositional inferences that can be 
made from the velocities and other data sug- 
gest that the crust is chemically different 
from the mantle and that it may have 
formed by the upward differentiation of the 
lighter minerals. 

Below the crust the lunar mantle may ex- 
tend to a depth of about 1000 km, probably 
with nearly constant compressional veloci- 
ties (refs. 81 and 82) . Data from moonquakes 
and meteorite impacts specify the properties 
of the mantle. The recording of well-defined 
shear waves from the deep moonquakes im- 
plies that the mantle must be sulficiently 
rigid to 700 or 1000 km depth to prevent 
appreciable shear-wave attenuation. 

Seismic shear waves from farside impacts 
and farside moonquakes that have pene- 
trated deeper into the Moon than about 700 
km are attenuated (refs. 82, 83, and 84). 
This indicates a possible “softening” of the 
material in the central region of the Moon 
below about 700 or 1000 km. Any softening 
that would reduce Q to less than about 500 
would explain the observations. This could 
be achieved by temperatures approaching 
the solidus, by a very small amount of partial 
melt, or by other mechanisms (such as per- 
haps an increase in the amount of volatiles 
in the deep lunar interior or a different bulk 
composition) . 

If we assume that the attenuation is due 
to temperature effects, this would place some 
limits, within the compositional constraints, 
on the temperatures below 700 km in the 
Moon. For a pyroxene-olivine-rich composi- 
tion this is 1600° C. If the Moon has an Fe- 
FeS-rich core (see ref. 85 for a discussion), 
the minimum temperature for it would be 
about 1000°C, given by the Fe-FeS eutectic. 
Thus, the present-day lunar “core” tempera- 
tures range from about 1000 to 1600°C. 

Constraints Implied by Lunar Viscosity, 
Seismicity, and Mascons 

The lunar interior described in the pre- 
vious section and shown in figure 5 can be 
divided into three units : a crust, a relatively 
rigid 1000-km-thick “lithosphere,” and a 
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relatively soft deep interior “asthenosphere.” 

The crust and the lithosphere have sup- 
ported lunar mass anomalies (“mascons”) 
for over 3 b.y. This requires that (1) at the 
time of the mascon formation and mare 
flooding the lunar lithosphere was thick 
enough to support the loads (ref. 86) and 
(2) the viscosity of the lithosphere must be 
greater than about 10^® poise (refs. 87 and 
88). This high viscosity requires tempera- 
tures significantly lower than the melting 
curve in the outer few hundred kilometers 
of the Moon. 

Another important constraint on the 
present-day thermal state is that the lunar 
interior is the very low level of seismic ac- 
tivity in the Moon. Total seismic energy re- 
lease in the Moon is about 10“ ergs/yr, 13 
orders of magnitude lower than that of the 
Earth (refs. 83, 89, and 90) . All moonquakes 
are very small, and those whose focal depths 
can be determined have hypocenters be- 
tween 700 and 1200 km depth. These indi- 
cate that if there is any convection or 
tectonic motion taking place in the Moon 
today, it must be below about 700 to 1000 
km. 

Chronology of Lunar Igneous Activity 

The time history of differentiation and ig- 
neous activity imposes the strongest con- 
straints on the early thermal state of the 
Moon. The chronology of the lunar igneous 
activity starts from the formation of the 
original crust about 4.6 b.y. ago as deduced 
from model ages of lunar rocks and soils 
(refs. 91 through 94). More recently a Rb- 
Sr age of 4.6 b.y. has been obtained for an 
Apollo 17 dunite fragment (72415) by Albee 
et al. (ref. 95) . To create a 60-km-thick feld- 
spathic crust described in the previous 
section requires the total or partial melting 
and differentiation of at least one-half of the 
total volume of the Moon (ref. 96). This 
requirement strongly constrains the initial 
temperature of the lunar interior. 

Following the formation of the initial 
crust, the chronology of lunar rocks (an ex- 



Fig^re 5 . — Schematic diagram of the structure of the 

Moon as inferred from lunar seismic data (ref. 

80 ). 

cellent summary has been compiled by Doe 
(ref. 97) ) is defined by three episodes in the 
lunar evolution: 

1. The period 4.4 to 4.0 b.y. ago was 
characterized by continuing igneous 
activity, brecciation by meteoritic 
bombardments, and metamorphism of 
the older rocks. Because of the subse- 
quent events, the detailed record of 
this episode has been largely obscured. 
Available data include ages of anortho- 
sites (4.1 b.y. for lunar rock 15415) 
(ref. 98), the formation ages of 
KREEP basalts (refs. 50 and 99), and 
a number of the “recrystallization” 
ages of the metamorphosed breccias 
from the Apollo 12, 14, 16, and 17 and 
Luna 20 sites (see ref. 97 for a com- 
pilation) . These are shown in figure 6. 

2. Episode of basin excavation (4.0 to 
3.9 b.y. ago) is characterized by in- 
tense bombardment, resetting of the 
geologic clocks, and excavation of the 
major basins such as Imbrium, Tran- 
quillitatis, Serenitatis, Nectaris, etc. 
(refs. 94 and 100). 

3. During the episode of mare flooding 
(3.9 to 3.16 b.y. ago), the basins 
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formed were filled with mare basalts. 
The ages of these lunar basalts are well 
documented by dating samples from 
the Apollo 11, 12, 15, and 17 and Luna 
16 missions (refs. 18, 92, 93, 101 
through 104, and others). 

The failure to discover lunar igneous rocks 
younger than about 3.16 b.y. is an extremely 
important constraint in the lunar thermal 
history. Since basalts were probably derived 
from sources to about 400 km depth, the 
termination of lunar magmatic activity rep- 
resents the end of melting or partial melt- 
ing in the upper mantle of the Moon. 

During the past 3.0 b.y., the lunar surface 
has continued to evolve as a result of mete- 
orite impacts and cratering. There have been 
local melting and metamorphism in this pe- 
riod. There is no direct evidence for extensive 
volcanic activity or basalt flooding in the 
past 3.0 b.y. 

Magnetism of Lunar Rocks 

Although the Moon does not have a mea- 
surable dipole magnetic field at present, the 
lunar crust is extensively magnetized. This 



Age , b. y. 


Figui’e 6 . — A summary of early lunar evolution and 
igneous activity at the surface, based on Rb-Sr, 
Ar^’-Ar'", and concordia ages of returned lunar 
samples. Samples from different missions are desig- 
nated by “A” for Apollo (A-11 is Apollo 11) and 
“L” for Luna. Sources for this summary are pre- 
sented in the text. 


remanent magnetization has been verified 
by magnetic field measurements from the or- 
biting Apollo subsatellites (ref. 105), on the 
surface (ref. 71), and from the remanent 
magnetization of returned lunar samples 
(refs. 106 through 109, and others). This 
stable remanent magnetization was most 
likely acquired on the Moon as the rocks 
cooled through the Curie point of iron in 
the presence of a magnetic field. The inten- 
sity of the ambient magnetic field at the time 
of the formation of these rocks is estimated 
to be greater than 1000 gammas, and prob- 
ably as high as a few thousand gammas 
(ref. 110). 

The origin of the ancient magnetic fields 
which magnetized the initial lunar crust and 
the lunar rocks is one of the most important 
and as yet unresolved problems in lunar stud- 
ies. From the ages of the crust and the rocks, 
the magnetizing field must have been pres- 
ent at least from 4.6 to 3.2 b.y. ago. The 
hypothesis of an external magnetizing field 
is not favored because of the difficulties as- 
sociated with having a steady field of 1000 
gammas or more near the Moon for a period 
of 1.5 b.y. 

The favored hypothesis of the internally 
generated magnetizing field would require 
either a lunar dynamo or magnetized lunar 
interior (refs. Ill, 112, and 113). In the 
dynamo model, it is necessary for the Moon 
to have a conducting iron or Fe/FeS core of 
sufficient size early in its history. The forma- 
tion of such a molten core requires high 
internal temperatures early in the lunar his- 
tory. In the case of a highly magnetized 
model of the lunar interior (refs. Ill and 
113), it is required that the temperatures 
below a few hundred kilometers depth re- 
main below the Curie temperature (about 
760°C), while all the magmatic activity and 
differentiation take place at shallower 
depths. Such a model would put severe re- 
quirements on the thermal evolution models 
and also require a high magnetizing field at 
the time of the lunar accretion. 

Thus, the magnetic history of the Moon 
according to the above models would require 
either extensive melting and differentiation 
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of the core at the onset of lunar evolution or 
temperatures less than about 800°C for the 
first 1.5 b.y. in the lunar mantle. 

TEMPERATURE MODELS 
FOR THE MOON 

In this section we present several specific 
models for the thermal evolution of the Moon 
which illustrate the effect of varying initial 
conditions and bulk radioactivities. Previous 
papers (refs. 23, 30, and 86) have discussed 
many of the models presented here and have 
assessed their ability to satisfy the con- 
straints listed in the previous section. In 
this section we divide the models into three 
subgroups : initially homogeneous, convec 
tive, and inhomogeneous models. 

Initially Homogeneous Models 

We begin by considering a rather unreal- 
istic model: initially cold (T = 0°C every- 
where). This is a valuable exercise in the 
respect that it provides a verification of the 
computational scheme. Figure 7 shows three 
present-day temperature profiles each for the 
given present-day bulk uranium concentra- 
tions. Clearly, none of the models shown 
here satisfies the constraint of near-surface 
melting in the first billion years of the 
Moon’s history. In fact, the basalt solidus 
(dashed line) is not reached at any point in 
time in the lunar evolution. Furthermore, 
while reasonable uranium concentrations of 
greater than 37 ppb provide present-day 
melting in the lunar interior, these models 
cannot generate enough heat to satisfy the 
above constraint. Since radioactive decay 
alone cannot provide the heat necessary for 
early melting and differentiation, we must 
consider models with high initial tempera- 
tures. 

In an earlier section it was stated that 
one might consider the effects of early heat- 
ing by tidal dissipation and/or solar wind 
flux as the Sun passed through a T-Tauri 
stage by investigating an initially com- 
pletely molten Moon. Such a model is shown 



Depth, km 

Figure 7. — Present-day temperature profiles for an 
initially cold Moon (0°C everywhere) , as a func- 
tion of present-day uranium concentration. The 
solidus of anhydrous mare basalt (ref. IH) is also 
shown. 

in figure 8. In this case, all the radioactive 
heat sources are immediately differentiated 
and concentrated toward the surface. The 
formation of a dense core (ref. 86) com- 
posed of iron or FeS would easily be ac- 
complished. Because of the depletion of heat 
sources in the lunar interior and the very 
effective heat transport in the convecting 
molten regions, this model cools rapidly and 
is completely solid after 2 b.y. At the present 
time it is completely solid and continuing to 
cool. With an average uranium composition 
of 60 ppb, the heat flow value is 29 erg/cm^-s, 
within the reported range of the Apollo 15 
and 17 measurements. Although this model 
satisfies the constraints of early initial melt- 
ing and heat flow, it does not allow for a 
small partially molten core as implied by the 
seismic data. 

Since the normal process of accretion can 
produce initial near-surface partial melting 
without employing some catastrophic event, 
we now consider a model of lunar thermal 
evolution with initial temperatures calcu- 
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Fignre 8 . — A model for thermal evolution as a func- 
tion of time in an initially molten Moon. Time, in 
billions of years since lunar origin, is indicated by 
the number adjacent to each temperature profile. 
On this and later figures, the Moon is partially or 
completely molten at those depths where the tem- 
perature profile lies along the solidus. 

lated according to equations (1) and (2). 
This then allows for early segregation of the 
lunar crust and differentiation of the radio- 
active heat sources. 

A model that satisfies all the major con- 
straints presently known is shown in figure 
9. The initial temperature is calculated from 
the accretion model, with an accretion time 
of 100 yr and a base temperature of 800°C, 
which allows for early heating by sources 
other than accretion. Initial melting extends 
to a depth of about 800 km, giving sufficient 
volume to differentiate a primordial lunar 
crust. The present-day uranium concentra- 
tion averages to 60 ppb. The thermal con- 
ductivity is taken to be constant at 0.45 X 
10® erg/cm-s-°C. 

It can be seen from figure 9 that this 
model cools at the surface, giving a solid 
lithosphere growing at a rate of about 220 
km per b.y. During the period of mare fill- 
ing, the upper extent of melting progresses 
from 140 to 280 km, in agreement with the 


depth of origin of mare basalts (ref. 114) 
and with the necessity of maintaining a 
thick, cool lithosphere to sustain the stresses 
associated with mascon gravity anomalies. 
Shortly after 1 b.y., a convecting core devel- 
ops due to undifferentiated heat sources of 
nearly primordial radioactive abundances 
near the center of the Moon. The core has a 
maximum radius of 1340 km at 1.5 b.y. and 
then slowly cools by conduction through the 
lithosphere. At present, this model is mol- 
ten below 1000 km. This is in agreement 
with S-wave attenuation in the core region 
and with the occurrence of nioonquakes at 
700 to 800 km depth. The surface heat flow 
is 30 erg/cm^-s, consistent with the mea- 
sured values. 

Convective Models 

We now consider an approximation to 
solid-state convection in the Moon. We pro- 
pose to simulate the efficient heat transport 



Figure 9. — Thermal evolution in a Moon accreted in 
100 yr at a base temperature of 800°C, with an 
average present-day uranium concentration of 60 
ppb. Symbols are explained in the caption to 
figure 8. 
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of convection by raising the thermal con- 
ductivity by a factor of 10 on a grid point 
that has a temperature greater than a critical 
temperature of 1000°C. This is equivalent 
to Tozer’s steady-state “core” tempera- 
ture. Otherwise, the conductivity follows 
the scheme of Schatz and Simmons (ref. 64). 
We permit this enhanced conduction after 
1 b.y., i.e., after rapid differentiation is com- 
pleted and a lithosphere has formed. As in 
the case of convection of molten material, 
this approximation does not include mass 
transport phenomena, nor is it a statement 
on any particular solid-state creep mecha- 
nism. It only enhances the heat transport. 

Figure 10 illustrates a model with convec- 
tion as described above. The initial condi- 
tions and present-day bulk uranium 
concentration are the same as for the pre- 
vious model (fig. 9). The evolution proceeds 
as before until 1 b.y. At this point, convec- 
tive processes dominate, and the conductivity 
is raised for temperatures above 1000°C. At 
the present, a convecting core of radius 1100 
km and temperature nearly 1000°C under- 



Fignre 10 . — The effect of convection by solid state 
creep on the thermal evolution of the Moon (see 
text). The initial parameters for this model are 
the same as for figure 8. 



Figure 11 . — Thermal evolution in a Moon derived 
from inhomogeneous accretion (ref. 30). The ini- 
tial temperature profile assumes a cold (0°C) ac- 
cretion over a 1000-yr time interval. Primary 
zoning of radioactivity is assumed (see text). 


lies a nonconvecting, cool, rigid lithosphere. 
The present heat flow is 29 erg/cm^-s. The 
features of this model are essentially the 
same as Tozer’s models. 

Inhomogeneous and “Blanketed” Models 

We have considered up until this point 
models which assume an initially homoge- 
neous distribution of heat sources in the 
Moon. The only exception was that of the 
molten Moon in which all the radioactive 
isotopes were initially differentiated to the 
surface. One type of model that allows 
relatively cool interior temperatures and sat- 
isfies the heat flow constraint without invok- 
ing solid-state convection has been suggested 
by several investigators (refs. 115, 116, and 
others). This involves a decrease in the con- 
centration of radioactivity in the Moon, with 
depth as a primary feature associated with 
inhomogeneous accretion. 

Figure 11 shoes a thermal history model 
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from Toksoz and Solomon (ref. 30) similar 
in many respects to the one proposed by 
McConnell and Cast (ref. 22). The initial 
temperature profile assumes accretion in 
1000 yr, with a base temperature of 0°C. We 
assume that the radioactivity is stratified 
such that the present-day uranium abun- 
dance is 120 ppb in the uppermost 200 km, 
50 ppb in the next 250 km, and 10 ppb in the 
remainder of the Moon. This gives a bulk 
uranium concentration of 53 ppb. The K/U 
and Th/U ratios are fixed with values given 
in table 1. The conductivity follows the 
Schatz and Simmons (ref. 64) model. Al- 
though the initial temperature barely ex- 
ceeds the solidus between 40 and 140 km 
depth, shorter accretion times would provide 
initial melting at greater depths as in the 
previous models presented in this paper. 
With differentiation allowed for in four dis- 
crete time steps between 0 and 0.9 b.y., the 
evolution of the upper few hundred kilome- 
ters is similar to that in figures 9 and 10. 
The present-day heat flow is 32 erg/cm^-s, 
and the interior of the Moon has tempera- 
tures well below the solidus. 

The difference between this model and the 
others discussed earlier is that the interior 
of the Moon is depleted in radioactive ele- 
ments at cold temperatures and, thus, heats 
slowly so that the Moon is everywhere solid 
at present. Whether or not a cold accretion 
can accompany primary chemical layering is 
debatable ; however, such an origin is incom- 
patible with a molten iron-rich core early in 
the Moon’s history. Furthermore, as in the 
case of solid-state convection, the present-day 
temperature profile does not seem to satisfy 
the constraint implied by the attenuation of 
S-waves in the deep interior. 

We finally examine the effect of a low 
thermal conductivity “blanket” layer on the 
Moon. Measurements of the effective thermal 
conductivity of uncompressed lunar soils 
(ref. 62) indicate values more than two or- 
ders of magnitude less than that assumed for 
the interior (see fig. 3). Some investigators 
have argued that the widespread distribution 
of these lunar fines following the large impact 
events 4 b.y. ago might have produced local- 


ized melting at shallow depths. Figure 12 
compares such a model with that of figure 9, 
100 m.y. after the soils are assumed to have 
been scattered.' Since we are interested in the 
maximum possible effect of such an event, we 
will assume the lunar regolith to extend to 
a depth of 10 km, with a thermal conductivity 
of 0.45 X 10^ erg/cm-s-°C. Both values are 
probably overestimates. The depth of dense 
fracturing in the lunar crust has been esti- 
mated to be less than 10 km by Toksoz et al. 
(ref. 80). Thermal conductivity increases 
rapidly with compression of the soils, and 
in-situ conductivity measurements made at 
the Apollo 15 and 17 heat flow sites (ref. 57) 
yield values only about 1 order of magnitude 
lower than the conductivity assumed for the 
bulk of the Moon. In any case, the maximum 
effect shown in figure 12 amounts to a rise 
in temperature of only about 180° C at 25 km. 
Thus, it can be confidently concluded that a 
uniform thermal “blanket” would have had 



Figfure 12 . — The maximum effect of a thermal 
“blanket” layer on near-surface lunar tempera/- 
tures 100 m.y. after the presumed time of the Im- 
brium event (600 m.y. after lunar origin). Curve 
B represents a model with initial conditions the 
same as figure 9, but with the introduction of a 
low-conductivity, 10-km layer at 600 m.y. Curve A 
is the near-surface temperature profile for the 
model in figure 9, with no blanket layer. 



312 


COSMOCHEMISTRY OF THE MOON AND PLANETS 


2000r^ — - 



Figure 13 . — Range of present-day temperature pro- 
files for theoretical thermal history calculations. 
The upper hound is from Toksoz and Solomon ( ref. 
SO), while the lower bound is the convection solu- 
tion of Turcotte et al. (ref. 66). 

a minimal effect on the early evolution of the 
Moon. 

While the models discussed above are lim- 
ited in number, they do serve to demonstrate 
the importance of initial conditions, geophys- 
ical constraints, and heat transport mecha- 
nisms in determining the thermal history of 
the Moon. 

The range of present-day temperature pro- 
files inside the Moon that satisfy the main 
constraints listed previously is shown in 
figure 13. One encouraging aspect of this fig- 
ure is that regardless of how the temperatures 
are calculated, whether conduction or convec- 
tion mechanisms are prevalent, the present- 
day temperatures are relatively narrow. 
Hopefully, with additional data, not only will 
the present-day temperature limits be nar- 
rowed, but also the early history of the Moon 
will be better understood. 

EVOLUTION OF THE MOON 

The thermal state of the lunar interior as 
a function of time is shown in figure 14, along 


with major episodes in lunar evolution. 
During the first 2 b.y., the lunar upper mantle 
undergoes sufficient melting to account for 
the differentiation of the crust and the sub- 
sequent lunar volcanism and mare filling. The 
zone of melting deepens with time. The litho- 
sphere thickens at the rate of about 220 
km/b.y. during this period. In particular, 
the shallowest melting progresses from 140 
to 280 km depth during the period of mare 
filling, in agreement with the depth of origin 
of mare basalts (ref. 114) and with the need 
for a reasonably thick lithosphere to sustain 
the stresses associated with mascon gravity 
anomalies. The disappearance of melting in 
the mantle coincides roughly with the termi- 
nation of magmatic events. At present the 
whole Moon is cooling. The deep interior, 
below a depth of 1000 km, may be hot enough 
for partial melting. This will explain the S- 
wave attenuation in this zone as described 
previously. 

If indeed partially molten, the rheologic 
properties of the deep lunar interior may be 

TIME (B.Y.) 
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Figure 14 . — Thermal evolution of the lunar interior 
as a function of time, based on the temperature 
calculations shotvn in figure 9. Major episodes of 
differentiation and igneous activity are shown on 
top and are based on lunar sample ages. Isotherms 
are in 0° C. Fine and coarse shading denote re- 
gions of partial and extensive melting, respectively. 
Zones of moonquake activity and high attenuation 
of S-waves shown on the figure apply to the 
present. 



EVOLUTION OF THE MOON AND TERRESTRIAL PLANETS 


313 


similar to the Earth’s asthenosphere and may 
have convective motions. Such convection 
may exert small stresses at the bottom of the 
lithosphere where the moonquakes occur, but 
it could not induce sufficient stresses in the 
thick lunar lithosphere to cause large moon- 
quakes or active tectonic motions. 

Present-day temperatures shown in figure 
13 and 14 generally exceed the melting tem- 
peratures of Fe or Fe/FeS combinations in 
the deep lunar interior. Thus, if there were 
a concentration of Fe or Fe/FeS in the cen- 
ter, it would be molten, and the Moon would 
have a molten core. 

Thermal Evolution of Mars, 
Venus, and Mercury 

In this section we calculate thermal history 
models for planets Mars, Venus, and Mercury. 
The Earth has been excluded from these cal- 
culations. Its evolution history is similar to 
that of Venus and is generally well under- 
stood. Its present-day temperatures are 
strongly constrained by the data, and models 
can be tailored to meet these requirements. 
New data that became available for Mars, 
Venus, and Mercury from recent planetary 
missions, however, justify a reexamination 
of the thermal histories of these planets. 

For the calculations we use the numerical 
techniques described previously and follow a 
procedure similar to that used for the Moon. 
Since there are many fewer data for the 
planets, the calculated thermal models are 
less rigidly constrained than the lunar mod- 
els. 

MARS 

The thermal history of Mars is not as 
strongly constrained as for the Moon or 
Earth. However, new data, primarily from 
the Martian Orbiter Mariner 9, provide in- 
formation that places some conditions on its 
evolution and internal structure. We will now 
present several theoretically calculated ther- 
mal evolution models from Johnston et al. 
(ref. 32). 


Structure and Constraints 

The two main pieces of data regarding 
the composition and structure of the interior 
of Mars are its mass (M) and moment of in- 
ertia factor {I /MR-). These roughly de- 
termine the density variation within the 
planet and may indicate the presence (or 
absence) of a high-density core. 

Using values of Jo and the mean radius as 
determined from Mariner 9 experiments 
(ref. 117), we calculate 7/MP- = 0.377 ± 
0.001, indicating a distinct density increase 
with depth. This increase could easily be ex- 
plained by a high-density core. Recent re- 
sults from the U.S.S.R. Mars probe indicate 
that Mars maintains a relatively weak dipole 
magnetic field (ref. 118), giving further cre- 
dence to the idea of a core. 

The presence of a core places a fairly 
strong constraint on the thermal evolution. 
For core formation to take place, tempera- 
tures within the planet must have been 
nearly everywhere above a core material 
melting curve at one time. If Mars evolved 
in a manner similar to that of the Earth and 
Moon, the differentiation of core and mantle 
most likely took place relatively shortly after 
origin. It thus seems to be the case with 
Mars that volcanism and, therefore, differen- 
tiation have been active since early in the 
Martian history (ref. 119). 

We assume that Mars has retained sulfur 
during its condensation and, thus, has an 
Fe-FeS core (refs. 58 and 120), and that 
potassium is partitioned into the sulphide 
phase (refs. 121 and 122), thus providing a 
heat source in the core by the decay of K^®. 
This model leads to early differentiation of 
Mars and early core formation, without re- 
quiring high initial temperatures necessary 
for the reduction of Fe or Si. This is possible 
due to the low and relatively pressure- 
independent Fe-FeS eutectic at 990°C (ref. 
123). 

Many new data concerning the evolution 
of Mars were obtained from the Mariner 9 
photographs. The general picture was a 
planet much more dynamic than was pre- 
viously thought. From the geologic interpre- 
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tation of these photographs (refs. 119, 124, 
125, and 126), generally four regions of 
prominence stand out: an ancient cratered 
terrain, younger volcanic plains, large vol- 
canic constructs such as shield volcanoes, 
domes, and craters, and extensive sedimen- 
tary deposits. 

The Martian evolution inferred from the 
above features can be summarized as follows 
(ref. 126) : 

1. Roughly 4.6 b.y. ago Mars accreted. An 
intense bombardment of the surface 
took place at, or somewhat later than, 
the final stages, resulting in the heavily 
cratered, oldest terrain. 

2. Volcanism, and probably the differen- 
tiation of the core, followed shortly 
thereafter. 

3. Intense volcanic and tectonic activity 
characterized by basalt plains, uplifted, 
faulted, and eventually eroded, was ap- 
parently the next stage of Martian evo- 
lution. It was in this stage that the 
volcanic shields such as Nix Olympica 
began to develop. 

4. The next stage featured large-scale 
tectonic activity and the completion of 
the shield volcanoes. 

5. The final stage of evolution, lasting 
until the present, shows evidence of 
moderate tectonic and volcanic activity 
and wind and possibly water erosion 
and deposition. 

Although no absolute ages are available for 
these stages of Martian evolution, they are 
valuable in constraining thermal history 
models. We know that an initial crust must 
have differentiated and the core formation 
must have occurred early in the Martian his- 
tory, with volcanism, implying mantle 
melting, beginning shortly after and proceed- 
ing to the present time. 

Heat Sources and Input Parameters 

The thermal evolution models of Johnston 
et al. (ref. 32) are calculated in the manner 
described previously, with the exception that 


the formation of an Fe-FeS core and the dif- 
ferentiation of K*® into the core are specifi- 
cally allowed. 

As before, the most important heat source 
in the thermal history calculations is the heat 
generated by the decay of the long half-life 
radioactive isotopes U“®, U^®, Th“*, and 
K'*®. No surface measurements of the abun- 
dances of these isotopes have been made for 
Mars, so we are forced to look at the only 
available data: those from the Earth, the 
Moon, and meterorites. The radioactive 
abundances for the Moon and some meteor- 
ites were discussed in an earlier section. The 
use of lunar, howarditic, or eucritic composi- 
tions in Mars, however, results in present- 
day temperatures too hot for any reasonable 
model. 

The “chondritic coincidence” for the Earth 
has long been recognized (ref. 3), giving 
U = 11 ppb, but this requires the depletion 
of potassium in the crust relative to chon- 
drites. Although putting the potassium in an 
Fe-FeS core can eliminate this problem, if 
all the excess potassium is put in the core, 
this increases the temperature significantly. 
In the other extreme, assuming no potassium 
in the core and K/U = 10 000, McDonald 
(ref. 7) recalculated the absolute uranium 
abundance for the Earth’s mantle, constrain- 
ing the heat production under steady-state 
conditions to match the average terrestrial 
heat flow. His uranium concentration, reduced 
to a homogeneous planet, is about 31 ppb. 
One might, therefore, place bounds on the 
Earth’s absolute abundances, with U from 
11 to 31 ppb and K/U from 80 000 to 10 000. 

Recent models of the chemistry involved in 
the condensation of the solar nebula (refs. 59 
and 127) indicate that a number of refrac- 
tory trace elements, including the rare 
earths, will crystallize along with Ca-Al-rich 
assemblages at about 1500 K. This is before 
the condensation of magnesium silicates and 
metallic iron. Potassium and other volatiles 
will condense at lower temperatures (between 
1100 K and 1200 K). Since these tempera- 
tures are higher than the condensation tem- 
perature for Mars (=« 450 K), we might 
expect that the planet retained solar propor- 
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tions of U, Th, and K. This amounts to a 
bulk uranium content of 15 ppb, with 
K/U = 50 000 and Th/U = 4. We have 
assumed, however, a bulk radioactive abun- 
dance equivalent to MacDonald’s Earth 
model; i.e., U = 31 ppb, and K/U = 10 000. 
It is important to note, however, that the 
total heat production of these abundances is 
the same as that obtained with solar propor- 
tions. The only difference in the thermal 
modeling arises in the percentage of potas- 
sium allowed to enter the sulphide phase of 
the core. For solar abundances, this neces- 
sarily requires less potassium than is assumed 
in this paper (to be discussed later), but 
since little is know quantitatively about the 
partitioning of potassium, and given uncer- 
tainties in solar abundances, we believe this 
to be a minor effect. 

The thermal conductivity of the Martian 
mantle is assumed to be temperature-depen- 
dent and described by the empirical equation 
derived by Schatz and Simmons (ref. 64). 
The silicate melting curve is taken to be an 
extrapolated peridotite solidus (ref. 128). 
Melting, convection of molten material, and 
differentiation of the radioactive isotopes are 
modeled as before. 

The formation of the Fe-FeS core begins 
when temperatures in some portion of Mars 
exceed the Fe-FeS eutectic of 990° C and is 
completed when the temperatures for the bulk 
of the planet are above the liquidus for an 
assumed composition of iron and sulfur. In 
the calculations, core formation is modeled 
as both a continuous and discrete event. When 
the temperatures exceed the eutectic at a grid 
point, a small amount of excess heat, equiva- 
lent to the mass fraction of core material, is 
diverted to the next lower grid point and con- 
verted to its temperature equivalent. Thus, 
the core material effectively melts its way 
into the interior. Also, a specified percentage 
of the available potassium is released from 
the silicates in one discrete step and allowed 
to enter the core. Internal heating is also 
slightly enhanced by the gravitational set- 
tling of the molten metal. The temperature 
field is reset to an adiabatic gradient when 
the internal temperature rises above the Fe- 


FeS liquidus where it is believed sufficient 
mass transfer has occurred to warrant this 
action. The gravitational separation of the 
core is irreversible, and, once formed, the 
core does not react with the mantle. 

We have assumed that the Fe-FeS eutectic 
and liquidus are unaffected by high-pressure 
phase changes of FeS. While this is appar- 
ently not the case (ref. 129), we maintain 
that the thermal models and conclusions pre- 
sented in this paper remain essentially un- 
changed. 

Parameters used in the thermal calcula- 
tions may be found in table 1. 

Thermal Models 

Figure 15, taken from Johnston et al. (ref. 
32), shows the effect of varying uranium 
concentrations. In the figure, the present-day 
temperatures in a uniform, initially cold 
{T — 0°C everywhere) Mars are shown as a 
function of the present-day concentration of 
uranium. The terrestrial value of K/U = 
10 000 is assumed. Also shown are the Fe- 
FeS eutectic, the Fe-FeS liquidus for 85 wt% 



Figure 15. — Present-day temperatures for an ini- 
tially cold ( 0° C) Mars, as a function of uranium 
concentration. Also shown are the Fe-FeS eutectic 
and liquidus and the extrapolated peridotite soli- 
dus used in this paper (ref. 32). 
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Fe, and the extrapolated peridotite solidus. 
No melting or differentiation is allowed. For 
U < 37 ppb, the temperatures are below the 
solidus. Clearly, these models do not satisfy 
the boundary condition that melting and dif- 
ferentiation of the silicates has occurred. For 
U = 20 ppb, the temperatures do not even 
exceed the Fe-FeS liquidus, prohibiting the 
completion of a core. To satisfy the boundary 
conditions, with the assumed radioactive heat 
production, higher initial temperatures must 
be used. 

A thermal model that satisfies the major 
constraints on the evolution of Mars is shown 
in figure 16 (model “B” of ref. 32). The ini- 
tial temperature is calculated from the accre- 
tion model, with an accretion time, t, of 10® 
yr and a base temperature of 500°C. The core 
begins to form immediately and is completely 
formed by 1.0 b.y. Dry silicates begin to melt 
at 2 b.y. This zone of partial melt expands 
and works its way into the interior. The sur- 
face cools slowly by conduction, so that at 
the present time the dry solidus is reached 
at a depth of 1800 km. The calculated heat 
flow for this model is 43 erg/cm*-s, which 
is between the average measured heat flows 
for the Earth and Moon. The core is pres- 
ently molten. 

The evolution of Mars and the extent of 
melting is shown schematically for this model 



Depth . ttm 

Figfure 16 . — Thermal evolution of a Mars accreted 
in 10‘ yr, vnth T*" = 500°C. U = SI ppb, with 
K/U = 10 000 (ref. 32). 
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Fi^re 17 . — The evolution of the Martian interior 
as a function of time, based on the temperature 
calculations shown in figure 16. Isotherms are in 
° C. The light shading above the dry silicate par- 
tial melt indicates the region where a small partial 
pressure of water would result in partial melt. 

in figure 17. Temperature and melting are 
plotted as functions of depth and time since 
origin. The light shading above the dry sili- 
cate melt shows regions where the introduc- 
tion of a small partial pressure of water 
would result in a partial melt. Two episodes 
of Martian evolution are apparent from the 
figure. The first is the early differentiation 
of a crust due to partial melting in the near- 
surface regions and the formation of an Fe- 
FeS core. The formation of the crust is 
crucial in the understanding of the surface 
geology. Early volcanism, evident from 
Mariner photographs and the age of many 
regions inferred from crater densities, im- 
plies the formation and maintenance of a 
crust early in the Martian history. This is 
followed by upper mantle melting and dif- 
ferentiation probably resulting in extensive 
surface volcanism. The region of mantle par- 
tial melt progresses toward the interior as the 
surface slowly cools by conduction. 

The presence of water and its relation to 
present-day tectonic activity is discussed in 
Johnston et al. (ref. 32). It was found that 
even for models such as the one discussed 
in this paper, where the anhydrous solidus 
is not reached until a depth of 1800 km, par- 
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tial melting can occur in the upper mantle. 
It is therefore concluded that one could ex- 
pect moderate tectonic activity at the present 
time. A 200-km-thick lithosphere obtained 
for this model is consistent with the depth 
required by hydrostatics to elevate magma 
to the summit of volcanoes like Nix 
Olympica — about 23 km above the mean 
Martian sphere — and with the gravitational 
roughness of the planet. 

Although the above temperature model is 
not a definitive statement about the present 
thermal state of Mars, it is adequate in the 
sense that all of the specified boundary con- 
ditions on thermal evolution are met. Hope- 
fully, data from the future missions will 
provide additional constraints, and, thus, 
thermal models will be better specified. 

VENUS 

Venus is perhaps the most unconstrained 
planet as far as theoretically calculated ther- 
mal models are concerned. It can be supposed, 
however, that from size considerations alone, 
the evolution of Venus might be similar to 
the Earth’s. There are, however, several im- 
portant differences. Condensation models of 
the solar nebula imply that Venus retained 
little or no sulfur, thus eliminating the pos- 
sibility of FeS or potassium in a core (ref. 
58). Therefore, proceeding on the bias of 
early differentiation, we must require high 
initial temperatures in order to form a core 
within the first billion or so years. If one 
takes the moment of inertia factor for Venus 
to be equivalent to the Earth’s (a rather 
strong assumption) it gives a core of radius 
about 2900 km for a composition of an Fe- 
Ni alloy. Lewis (ref. 58) also suggests that 
this core is surrounded by a massive mantle 
of Fe”-free magnesium silicates. A crust, 
similar in composition to the Earth’s, might 
also be expected. 

Radar maps of the Venusian surface indi- 
cate relief of several kilometers (ref. 130). 
Considering the density and velocity of the 
atmosphere, we believe these features to be 
subject to extensive erosion and therefore 


young and probably of tectonic origin. David- 
son and Anderson (ref. 131) have proposed 
that the rate of volcanism and tectonic activ- 
ity is greater on Venus than on the Earth. 
Weertman (ref. 132), however, believes this 
not to be the Case. Resolution of this conflict 
might be found from thermal models. Cer- 
tainly, though, if the stresses associated with 
the nonhydrostatic state of Venus are not 
supported statically, we would have to expect 
convective motions in the upper mantle. 

We now present a theoretical thermal 
model for Venus (fig. 18) . The computational 
parameters for this model may be found in 
table 1. The initial temperature profile is cal- 
culated from equations (1) and (2), with an 
accretion time of 0.25 X 10® yr taken to en- 
sure core differentiation within 1.0 b.y. The 
effect of adiabatic compression on the initial 
temperature is also included. The total heat 
production due to radioactive decay is essen- 
tially chondritic, having assumed urahium 
and potassium abundances equivalent to 
MacDonald’s (ref. 7) Earth mantle values, 
i.e., U = 30 ppb, and K/U = 10 000 at the 
present day. This is consistent with theories 
on the condensation of the solar nebula (ref. 
59). The bulk uranium concentration for 



Figure 18 . — A model for the thermal evolution of 
Venus. The accretion time is 0.25 X 10' yr, and 
the effect of compressional heating is included in 
the initial temperature. Input parameters for this 
model are described in the text. 
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Venus might actually be slightly higher than 
the Earth’s, due to the fact that it has not 
retained sulfur and, thus, has a lower mean 
atomic weight (ref. 58). The thermal con- 
ductivity of the undifferentiated planet and 
of the differentiated mantle is from the 
Schatz and Simmons (ref. 64) model. 

Core formation is initiated when the tem- 
peratures within the interior exceed the Fe 
melting curve calculated for a homogeneous 
planet (ref. 133). At this point (about 1.0 
b.y.), it is assumed that the rate of inflow 
of molten iron to the core is sufficient so that 
the gravitational energy due to core separa- 
tion (refs. 134 and 135) allows for immediate 
core formation and temperatures in the in- 
terior are raised by the heat equivalent. We 
reset the temperature profile to an adiabatic 
gradient after core formation, due to the 
extensive mass transfer that has presumably 
taken place. Radioactive isotopes originally 
in the core region are displaced and distrib- 
uted homogeneously throughout the mantle. 

The temperature rise due to core infall is 
sufficient to initiate mantle melting immedi- 
ately after core formation. The mantle 
melting curve is derived from Lindemann’s 
equation which can be written in a form so 
that T„ = const where is the melting 
temperature and is the seismic parameter 
(D. H. Chung, personal communication). 
Using measured values of <f> for the Earth, 
this equation provides an adequate estimate 
for the mantle melting curve for Venus. Con- 
vection of molten material and differentia- 
tion of the radioactive isotopes are modeled 
as described previously. 

Present-day temperatures are shown in the 
curve labeled 4.6 b.y. The iron melting curve 
shown in figure 18 is calculated for a differ- 
entiated planet. Notable features of this 
model are the thin (100-km) lithosphere, 
partially molten upper mantle, solid lower 
mantle, and molten core. Despite the molten 
core, the absence of a magnetic field for 
Venus can be explained by the lack of heat 
sources in the core and the slow rotation rate 
of the planet. Calculated heat flux is 98 
ergs/cm^-s. 

Figure 19 shows schematically the evolu- 


tion of Venus as a function of time. Isotherms 
are in degrees centigrade. As noted before, 
core formation takes place at 1.0 b.y., and 
present-day conditions indicate a region of 
high fractional partial melt (shaded area) 
in the upper mantle (higher than expected 
in the Earth) , overlain by a lithosphere about 
100 km thick. Although the degree of partial 
melt may be a function of the high surface 
temperature (presumably due to the atmo- 
sphere which may not have been in existence 
throughout the history of Venus), we have 
computed thermal models with the surface 
temperature arbitrarily set at 0°C. This re- 
sults in a change of partial melt only in the 
upper few grid points. We therefore believe 
that any stresses associated with hydrostatic 
disequilibrium are supported within Venus 
by strong convective motions in the upper 
mantle. This points directly to the possibility 
of plate motion, similar to the Earth’s and 
likely at a higher rate, and leads us to expect 
significant tectonic activity at the present 
time. 

MERCURY 

The knowledge of the planet Mercury has 
increased significantly since the Mariner 10 
mission of April 1974. Preliminary results 
of television reconnaissance of the planet’s 



Figure 19 . — The evolution of Venus as a function of 
time, showing core formation and mantle melting 
( shaded region). Isotherms are in ° C. 
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surface (ref. 136) show it to be remarkably 
similar to the lunar highlands. Regions re- 
sembling lunar mares have been observed, 
implying volcanic activity during Mercury’s 
history. The presence of a magnetic field (ref. 
137) could imply a molten iron core. We pro- 
pose, then, that Mercury, as the other ter- 
restrial planets, differentiated early, forming 
a crust similar to the Moon’s and an iron 
core. 

Siegfried and Solomon (ref. 33), using 
methods similar to those described earlier in 
this paper, have calculated several thermal 
and density models for Mercury. They point 
out, however, that Mercury is different from 
the other planets in the respect that conden- 
sation models of the solar nebula imply en- 
richment in refractory elements such as 
uranium and thorium and retention of little 
or no potassium and other volatiles. Further- 
more, the thermal evolution is dominated by 
the high conductivity of iron, which makes 
up approximately 69 percent of the planet by 
weight. 

Using appropriate initial temperatures, 
Siegfried and Solomon have shown that 
Mercury can indeed differentiate and form 
a core. Difficulties arise, however, if one as- 
sumes no heat source in the core and differ- 
entiation takes place earlier than about 3.0 
b.y. These difficulties occur because the core 
is presently solid, with temperatures several 
hundred degrees lower than the iron melting 
curve calculated following Higgins and Ken- 
nedy (ref. 133). This is in apparent conflict 
with the presence of a magnetic field, if such 
a field is indeed due to an active core dynamo. 
Furthermore, differentiation at times earlier 
than 1.0 to 1.5 b.y. seems to be required from 
the nature of Mercury’s surface. 

We have therefore calculated a thermal 
model for Mercury (fig. 20), assuming the 
presence of some sort of heat source in the 
core. We have chosen this to be, as a matter 
of convenience, the decay of K^». The initial 
temperature profile is calculated from the ac- 
cretion model, with t = lO^ yr, and T^, = 
1127°C (1400K), which provides initial iron 
melting to a depth of about 900 km. Thus, 
core formation and mantle melting are ini- 


tiated immediately. Other model parameters 
such as uranium and thorium abundances are 
taken from Seigfried and Solomon (ref. 33) 
and may be found in table 1. Mantle conduc- 
tivity is taken to be constant and equivalent 
to that used for the lunar thermal model 
presented in this paper {K = 0.45 X 10« erg/ 
cm-s °C). 

As shown in figure 20, temperature profiles 
indicate a solid lithosphere 200 km thick at 
2 b.y. and about 500 km thick at the present 
time. The core is presently partially molten 
and therefore might be able to sustain con- 
vective motions required to maintain a mag- 
netic field. The lithosphere is probably thick 
enough to ensure a low level of tectonic ac- 
tivity. We would conclude from this model, 
however, that the regions on Mercury’s sur- 
face that appear like the lunar mare may 
actually be younger, with tectonic activity 
having been prevalent on the planet for a 
longer period of time. We have also found 
that a concentration of 0.156 ppm at the 
present time provides a minimum value re- 
quired to maintain a partially molten core. 
The calculated heat flux is 60 ergs/cm®-s. 

For comparison we show a thermal model 
of Mercury in figure 21 that contains no heat 
source in the core. Other model parameters 
are the same as those described above. In 
this case, although the mantle and core dif- 
ferentiate, the entire planet has temperatures 
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Fignire 20 . — A thermal model of Mercury, with a 
heat source assumed in the core. Input parameters 
for this model are described in the text. 
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origin Time , b.y. present 

Figure 21 . — The thermal history of Mercury, with 
the same parameters as for the model in figure 
20, except that the core is assumed to be free of 
heat sources. 
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Figure 22 . — The evolution of Mercury as a function 
of time. Shading above the core region indicates 
mantle melting. Isotherms are in " C. 

below the melting curve at the present time. 

Figure 22 shows the evolution of Mercury 
as a function of time for the model in figure 
20. Given the parameters used in this model, 
we conclude that Mercury could have a mol- 
ten iron core if internal heat sources are 
assumed, and that present-day temperatures 
and the thickness of the lithosphere limit the 
tectonic activity to very low levels, if there 
is, in fact, any activitv. 

Discussion and Conclusions 

The thermal history models calculated for 
the Moon, Mars, Venus, and Mercury are 


representative of the whole spectrum of sizes 
and compositions of the terrestrial planets. 
The evolution histories are shown together in 
figure 23 for a comparative analysis. It is 
clear from this figure that although the ini- 
tial temperatures are similar for these plan- 
ets the evolution histories and present-day 
temperatures are quite different. The size 
(radius) of the planet more than any other 
factor seems to control its evolution and 
present-day thermal state. For the Moon and 
Mercury, the two smallest bodies, the evolu- 
tion and differentiation took place relatively 
early in their histories. At present, both have 
thick lithospheres and are mostly solid. 
Venus, on the other hand, seems to be still 
at the peak of its evolution. Mars falls in 
between. 

A good illustration of the comparison of 
evolutionary history of the planets can be 
obtained from the value of the total thermal 
energy of the planet as a function of time 
given by: 

E{t) = iUPC^T(t)dV + Q{t))/M 

where E(t) is energy in ergs per gram, T{t) 
the absolute temperature, and Q(t) the heat 
that goes into fusion at time t. M is the mass, 
and V the volume of the planet. This is il- 
lustrated for the planets examined in this 
paper in figure 24. 
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Figure 23 . — A comparison of the evolution of the 
planets. 
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Figure 24 . — Total thermal energy as a function of 
time for the Moon, Mars, Mercury, and Venus. 


Several conclusions may be drawn from 
this figure. First, the maximum energy cor- 
responds to the peak of the evolutionary pro- 
cess. Second, the rate of cooling is dependent 
primarily on the planet’s size and somewhat 
on composition (that is, thermal conductiv- 
ity). 

For example, the Moon, which is the small- 
est planet, has a peak of thermal energy 
during the first billion years after formation. 
This corresponds roughly to the period of 
large-scale magmatic activity as evidenced 
by the mare basin filling. Afterwards, the 
energy level decreases rapidly as the Moon 
cools. At the present time we would expect 
the Moon to be relatively inactive tectoni- 
cally. This is illustrated by a simple compari- 
son with the Earth, which appears to be at 
the peak of its evolution at present. The ter- 
restrial lithosphere is about 100 km thick, 
while that of the Moon is about 700 km. The 
Earth is tectonically active, with seismic en- 
ergy release of about E, = 10*® ergs/yr, 
while the Moon is inactive with E, ^ 10^^ 
ergs/yr (ref. 90). While three-fourths of the 
Earth’s surface is covered with rocks of 500 
m.y. or younger, the youngest crystalline 
rock found on the Moon is 3.16 b.y. 

Mercury, the next smallest planet, also 
reaches its peak of evolution early in its his- 


tory. We would expect, upon examining figure 
24 (with thermal energy calculated for the 
model shown in figure 20), that the bulk of 
the tectonic activity occurred before 2.0 b.y. 
The sharp increase in thermal energy is due 
to core formation. After 2.0 b.y. the energy 
decreases, showing that the planet is pres- 
ently cooling and is probably dormant. If 
one does not include heat sources in the core, 
as in figure 21, the slope of E(t) would be 
greater at the present time and the planet 
would obviously be cooler, heat transfer 
being dominated by the high conductivity of 
iron. 

Mars, intermediate in size among the ter- 
restrial planets, reached its peak of evolution 
around 3 b.y., maintained a relatively con- 
stant level of thermal energy for about a bil- 
lion years, and is now beginning to cool. From 
this we can conclude that, while magmatic 
and tectonic activity have existed on Mars 
since early in its history, most of the volca- 
nism and tectonics probably occurred in the 
last 2 b.y. We also expect a moderate level 
of tectonic activity at the present time. 

Venus, in many respects similar to the 
Earth, differentiated a core early and has 
maintained a fairly constant level of thermal 
energy throughout its history. It can be con- 
sidered to still be at its peak in evolution. 
Thus, we can expect tectonic processes in 
Venus similar to those observed on the Earth. 

In interpreting the above results, one must 
be careful to note that the thennal energies 
calculated for each planet were taken from 
one nonunique thermal history model. These 
models, however, are reasonable in that they 
satisfy most of the constraints on the evolu- 
tion of the planets known at the present 
time. Since they represent a bias toward early 
differentiation of all the terrestrial planets, 
the thermal energy turns out to be a useful 
tool to examine the effect of planetary size 
on the durations of high evolutionary and 
therefore tectonic activity. 

In summary, we conclude the following: 

1. The terrestrial planets were heated in 
some manner during or shortly after 
accretion. 
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2. This initial heating allowed for early 
core formation and/or differentiation. 
Early differentiation is evident from 
rock ages for the Earth and Moon and 
from photographic interpretation of 
Mars and Mercury surface features. 

3. For planets with large cores such as 
the Earth, Venus, and Mercury, the 
energy released upon core separation is 
a large portion of the total heat budget 
of the planet. 

4. The Moon is characterized as a dif- 
ferentiated body which evolved early 
in its history. The models that satisfy 
the constraints define a relatively nar- 
row range for the temperatures in the 
outer several hundred kilometers, re- 
gardless of conduction and convective 
type calculations. Below 1000 km, the 
control on temperatures is not as great. 
Temperatures at the center of the Moon 
may be anywhere between 1000° and 
1600° C. At these temperatures the 
Moon could have a partially molten (or, 
if Fe and FeS are differentiated, a to- 
tally molten) core. 

Acceptable evolution models require 
extensive differentiation and upward 
concentration of radioactive heat 
sources early in lunar history. This dif- 
ferentiation most likely is a product of 
extensive melting, although differentia- 
tion due to inhomogeneous accretion 
cannot be ruled out. A model of the 
Moon which is initially hot and exten- 
sively or totally molten can best satisfy 
the constraints. 

5. Thermal models for Mars, assuming an 
Fe-FeS core composition, indicate for- 
mation of the core and the differentia- 
tion of an early crust within the first 
billion years after origin in order to 
satisfy the given constraints on the 
thermal evolution. These models in- 
volve large-scale differentiation of the 
mantle silicates in the past 2 b.y. This 
would be accompanied by volcanism 
and outgassing of volatiles at the sur- 
face. The inclusion of water in the up- 
per mantle system leads to partial 


melting at a depth of about 22 km or 
greater at the present time. 

6. Venus, while not as strongly con- 
strained as the Moon or Mars, is char- 
acterized as a planet not unlike the 
Earth in many respects. Core forma- 
tion is allowed to occur by 1.0 b.y., and 
present-day temperature profiles indi- 
cate a partially molten upper mantle 
overlain by a lithosphere approxi- 
mately 100 km thick and a molten 
Fe-Ni core. 

7. Mercury is taken to be similar to the 
Moon except for its large iron core. 
Thermal models with an assumed heat 
source in the core yield a planet that 
has a solid mantle to about 500 km 
depth, with a partially molten core. 
Those with no heat sources in the core 
have present-day temperatures well 
below the melting curve for iron, heat 
transfer being dominated by the high 
thermal conductivity of the iron. 

8. The examination of total thermal 
energy as a function of time for each 
planet provides a method for determin- 
ing the effect of planetary size on the 
durations of high tectonic activity. 
The Moon, smallest in size, is presently 
cool and tectonically inactive. Mercury 
is probably inactive at the present 
time, but its peak of activity lasted 
longer than the Moon’s. Mars, interme- 
diate in size, reached its peak in evolu- 
tion in the past several billion years, 
but is probably moderately active at 
the present time. Venus, having been 
able to retain its thermal energy 
throughout its history, may today have 
tectonic processes similar to the 
Earth’s. 
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